Optimizing the efficiency of evaporative cooling in optical dipole traps 
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We present a combined computational and experimental study to optimize the efficiency of evap- 
orative cooling of optically-trapped atoms. We employ a kinetic model of evaporation to determine 
the optimal relation of trap depth to average trap frequency given experimental initial conditions. 
We then experimentally implement a highly efficient evaporation process in a misaligned crossed- 
beam far off-resonance optical dipole trap (MACRO-FORT), showing excellent agreement between 
the theory and experiment. This method has allowed the creation of Bose-Einstein condensates of 
87 Rb with 2xl0 4 atoms starting from only a few 10 5 atoms initially loaded in the optical trap. 

PACS numbers: 64.70.fm, 67.85.Hj 



Over the last 30 years, the control and study of ul- 
tracold atoms has continued to impress. Initial studies 
in both optical [1] and magnetic [2] confinement of neu- 
tral atoms led to the realization that evaporative cooling 
could create colder atom temperatures [3 . Subsequently, 
evaporative cooling played a critical role in realizing 
atomic Bose-Einstein condensates (BECs) [IHS] and de- 
generate Fermi gases [7] . Though primarily used as a tool 
in the creation of quantum degenerate gases, evaporative 
cooling also presents interesting physics in itself [5HI3]. 
The goal of most evaporative cooling schemes is to opti- 
mize the number of atoms, N, that reach quantum degen- 
eracy. Forced evaporation proceeds by lowering the trap 
depth, [/, allowing the hotter atoms to escape the trap 
while the colder atoms remain and rethermalize. The ef- 
ficiency of evaporation is characterized by the parameter 
leff = -log(pf/pi)/log(N f /Ni), where p = n\ 3 dB is the 
phase space density, n = A/^ 3 [m/(27r/c#T)] 3 / 2 is the cal- 
culated peak atomic density, X^b — \/2irh 2 /mkBT is the 
thermal deBroglie wavelength, T is the atoms' temper- 
ature, uj is the average trap frequency, m is the atomic 
mass, ks is Boltzmann's constant, and H is the reduced 
Planck's constant. 

Atom losses via three-body combination and one-body 
background collisions limit the attainable efficiency in 
most experimental evaporative cooling processes [T4] . 
The rates of these losses, T^b and Tib respectively, 
should be kept small compared to the rate of evapora- 
tion, T ev . The evaporation rate depends on the elas- 
tic collision rate of the atoms, T e i, and the ratio of the 
trap depth to atom cloud temperature r] = U/ksT. Us- 
ing a truncated Boltzmann distribution approximation, 
T ev = XT e i where A « (77 — 4)e _?7 for a 3D harmonic 
trap when rj > 6 [9J fl5] . The elastic collision rate of 
the atoms in a 3D harmonic trap is T e i = nav/ (2\/2), 
where a = Sna 2 is the elastic cross section, a s is the 
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s-wave scattering length (a s = 98ao for 87 Rb where ao 
is Bohr radius [II]), and v = ^{ksT /irm) 1 ^ 2 is the av- 
erage relative velocity of the atoms. The three-body 
re-combination rate depends on the density distribu- 
tion and atomic species: T^b = L J n(f) 3 df/N , where 
L = 4.3(±1.8) x l(T 29 cm 6 /s for 87 Rb in the F = 1 ground 
state [T7j . 

Motivated by the accessibility of all magnetic spin 
states, less stringent vacuum requirements, and the abil- 
ity to use Feshbach resonances to tune atomic interac- 
tions, all-optical trapping is increasingly popular in cold 
atoms experiments. The simplest and most commonly 
used evaporative cooling in an optical trap is done by 
simply lowering the trap laser power, and thus the trap 
depth. For such evaporation, the average trap frequency 
is reduced as uj oc U - 5 [TEJH9]. Since the collision rate 
(and thus evaporation rate) decreases with decreasing 
a;, efficient all-optical evaporative cooling can be diffi- 
cult due to stagnation of the cooling process. To address 
this limitation, various techniques have been developed 
to change the relationship of the trap frequency and trap 
depth (uj (x U u ) during the evaporation, such as the zoom 
lens trap [20] , tilted trap [I5j [21] , and adding additional 
lasers for confinement [22] . While these techniques aided 
in the achievement of BECs, reported values of j e ff have 
been between 2.5 and 3.5. Little attention has been paid 
to optimizing the efficiency of the evaporative cooling by 
optimizing v. In this work, we optimize the efficiency of 
evaporation by optimizing both r\ and v and experimen- 
tally demonstrate a j e ff = 4.2. 

We achieve high efficiency evaporation by choosing an 
evaporation route that maximizes the evaporation rate 
while keeping three-body losses low. We experimentally 
implement this route using a misaligned crossed beam 
optical dipole force trap (MACRO-FORT), which allows 
for tailoring v in the uj oc U v relation [2T] , In order to 
find the optimal route our experimental setup, we model 
the evaporation by numerically solving the energy and 
atom number of the atomic cloud as functions of time. 
For our model, we follow the kinetic theory of Ref. [9] 



which assumes sufficient ergodicity in the trap (see also 
Refs. [I3j [I5j |23] , noting [24]). We also assume operation 
in the deep trap (77 > 6) regime such that the density dis- 
tribution is approximated by n(r) = nexp[— (/(rj/ksT], 
and trap potential by U (r) = \m (oo 2 x 2 + u 2 y 2 + uj 2 z 2 ) . 
With this approach, E — 3NksT and the energy and 
atom number evolution equations take the form: 



E = -NT ev (r] + n)k B T + vE- r 



T 1B E-T 3B -E (1) 



N : 



-(r ev + r 1B + r 3B )N 



(2) 



where k « (77 — 5) / (77 — 4) for a 3D harmonic trap. In 
Eqn. [I] the first term accounts for the energy removed 
by evaporative cooling, the second term for the adiabatic 
weakening of the trap [25], the last two terms for the 
energy lost by one-body (e.g. background gas collisions) 
and three-body processes, respectively. Using this model 
and given the initial conditions after loading atoms in 
the optical dipole trap, the ideal values of 77 and v can 
be found to maximize J e ff- 
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FIG. 1. (a) Numerical solution of the efficiency of evapora- 
tion, Jeff , as a function of 77 and v given our experimental 
initial conditions. The solution is presented for parameters 
of 77 and v to achieve pf = 0.215 with a total evaporation 
time T evap < 10 s (area bounded by the thick, dotted line). 
The maximum efficiency, Jeff,maxi 

for these conditions is 4.51 
(square). Our experimental implementation employed values 
of 77 and v that achieved j e ff = 4.2 (circle). Initial conditions 
of evaporation are T = 60 pK, N = 5 x 10 5 , uj = 2tt x 1000 
Hz, Tib = 1/12 sec -1 . Since Ni, pi and pf are fixed, j e ff 
depends only on Nf. (b) and (c): Effect of varying the initial 
atom number, Ni, while keeping other initial conditions the 
same as in (a) on the optimal value of v for evaporation, (b), 
and resulting Jeff, max, (c). The grey dashed line marks the 
value of Ni in our experiment. 

Fig. [IJi shows the calculated j e f f given our initial con- 
ditions and evaporating to near the critical temperature 
T c (defined from k B T c = ^(iV/1.202) 1 / 3 ) for different 
values of 77 and v. We find a region of optimal 77 and nu 
centered around 77 = 8.95, v = 0.05. For our MACRO- 
FORT setup, limitations in the laser power available re- 
strict v to be between 0.15 and 0.50. The large tolerance 
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FIG. 2. Calculated trapping potentials for various MACRO- 
FORT powers along three axis (a-c) and a schematic of the 
experimental setup (d). The calculated trapping potentials 
in red, blue and black are shown in (a-c) for various com- 
binations of the optical trap beams' powers, corresponding 
tO (P wlde ,Pnarro W ) =(12W, 0.09W), (3W, 0.058W), (1.25W, 
0.20W) respectively. The 45 pm offset between the wide and 
narrow beams in the y-direction, marked by two-dotted lines, 
can be seen in (b). For (d), the magneto optical trap mag- 
netic field coils are in green, dipole trapping beams in orange, 
MOT trapping beams in red, and imaging beam in purple. 
Not pictured are the bias coils, which are used to cancel the 
earth's magnetic field and apply constant offset fields. Grav- 
ity is along the positive y-axis. 



of the parameters when determining j e f f allows for some 
flexibility in the experimental implementation of 77 and 
v, and we have designed our experimental evaporation 
ramp to have 77 ~ 8.6 and v = 0.22, so as to prefer- 
entially avoid three-body losses while still maintaining a 
high Jeff- As seen in Fig. [lb and c, loading more or less 
atoms initially in the dipole trap changes the ideal v as 
well as the maximum j e f f (while 77 is minimally affected 
by Ni). We note that we have assumed constant 77 and 
v for optimizing the evaporation process. At the cost of 
additional complexity, 77 and v can be varied during evap- 
oration to achieve slightly more efficient evaporation. For 
the initial conditions of our experiment, we find that im- 
plementing a variable v(t) during evaporation gives only 
a small gain (^0.1) in the numerically computed J e ff- 

To experimentally generate the MACRO-FORT, we 
use a single-frequency, single spatial mode, Erbium fiber 
laser with wavelength at 1550nm (IPG Photonics ELR- 
50-1550-LP-SF). The wide dipole trap beam has a beam 
waist of 88 pm and 18 W of initial power (P W ide)- The 
narrow beam (used to produce the "dimple" potential) 
has a waist of 20/im, 9W of initial power (P na rrow), 
crosses the wide beam at an angle of 66°, and is offset 
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45 /im radially from the wide beam's focus (see Fig. [2]). 
To load the MACRO-FORT, we follow a procedure simi- 
lar to Ref. [21] , which was developed to allow for the load- 
ing of a dipole trap from a magneto-optical trap (MOT) 
even when the dipole trap laser creates a strong AC Stark 
shifts on the cooling transition of the MOT (up to ~ 150 
MHz blue-shifts for our experiment) [26]. When loading 
the MACRO-FORT [27 , if only the wide beam is ap- 
plied we observe n = 1 x 10 11 atoms/cm 3 and T = 10 fiK. 
For such initial conditions, the T e \ is too low for efficient 
evaporation. The addition of the dimple potential after 
loading increase T e i , aiding efficient evaporation [28U30] . 
After 500 ms hold time in the MACRO-FORT with the 
dimple on, there are 5 x 10 5 atoms at 60 jaK and a density 
of 3.5 x 10 12 atoms/cm 3 . 
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FIG. 3. (a) Measured average trap frequency (Co) verses trap 
depth (U) for the near-optimal evaporation trajectory used in 
our experiment. The trap frequencies were measured by the 
parametric heating method. The experiment used v — 0.22 
for the first 3 seconds of evaporation. For comparison, the 
dashed line indicates the uj oc U° ' 5 relationship used in con- 
ventional optical evaporative cooling, (b) Extracted 77 during 
evaporation based on the T of atoms measured from time-of- 
flight, and trap depth U calculated (assuming ideal gaussian 
beams) from the measured dipole trap laser powers and beam 
geometry. 

We show the actual 77 and v during the evaporation 
ramp in Fig. [3] Forced evaporation with v — 0.22 and 
77 « 8 — 10 (noting the relative insensitivity of j e ff to 
small deviations of 77 and v from the optimal target val- 
ues) was implemented by exponentially decreasing the 
power in both beams over the course of 3 seconds to 12 W 
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FIG. 4. The experimentally measured atom number (a), 
atomic density (b), and phase-space density (c), as functions 
of time during evaporation. The theoretical calculations uti- 
lize the model of Eqns. [I] and [2] for v — 0.22, 77 = 8.6, 
uJi = 2tt x 1000 Hz, and assuming our experimental initial 
conditions. The slight overestimate of the theoretical n and 
p during the first 2.5 seconds is likely due to the presence 
of atoms in the "wings" of the narrow beam trap (e.g. see 
Ref. [22]). The insets in (c) show resonant absorption images 
(taken after 10 ms of time-of- flight expansion) of the atomic 
cloud at three different times along evaporation ramp, show- 
ing the thermal cloud to BEC transition. 



and 0.068 W respectively, achieving a high efficiency of 
j e ff = 4.2. At this point, there are 10 5 atoms at just over 
1 11K (see Fig.0. The average trap frequency is 475 Hz 
so Tq ~ 900 nK. Proceeding to evaporate in such a tight 
trap, however, would yield high 3-body losses. Thus, over 
the final 1.35 seconds, the power in the wide beam is low- 
ered to 1.25 W and the narrow beam to 0.020 W. This 
results in a nearly pure BEC of 2 x 10 4 atoms in a trap 
with measured frequencies of 27rx (60, 100, 270) Hz. Dur- 
ing this final 1.35 seconds, due to the atom cloud being 
near the BEC transition, the evaporative process is aided 
by the bosonic enhancement [3TJ [32] and no longer well 
described by the classical evaporative theory of Eqns. [l] 
and [2] This yields efficiencies even higher than 4.2. Af- 
ter further experimental optimizations, we increased the 
initial number of atoms loaded in the dipole trap and 
achieved pure BECs of 4 x 10 4 atoms. 

In summary, we have described a scheme to optimize 
the efficiency of evaporative cooling, where we have par- 
ticularly tailored the relation of the trap frequency to 
trap depth during the evaporation process. As a result, 
we achieve a j e ff of 4.2, the highest reported value to 
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our knowledge for all-optical evaporation. Different from 
Refs. [I5j [I9j l2T] . we include atom losses in modeling 
the evaporative process and optimize both 77 and v for 
our experimental conditions. Different from the evapo- 
ration optimization scheme in Ref. [13 , we optimize the 
evaporation by tuning v for the whole evaporation ramp, 
rather than changing 77 during evaporation. In addition, 
we have highlighted the utility of MACRO-FORTs that 
allow selecting the 77 and v values for all-optical evapo- 
ration and note that BECs can be obtained even when 
starting from a small number (a few 10 5 ) trapped atoms. 



The method shown here of optimizing the efficiency by 
finding optimal 77 and v is not specific to 87 Rb, and can 
be applied to achieve high j e ff for any atomic species 
that is evaporatively cooled in optical traps. 
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